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Recent experiments on similarly shaped polymer micro-cavity lasers show a dramatic difference 
in the far-field emission patterns. We show for different deformations of the ellipse, quadrupole 
and hexadecapole that the large differences in the far-field emission patterns is explained by the 
differing ray dynamics corresponding to each shape. Analyzing the differences in the appropriate 
phase space for ray motion, it is shown that the differing geometries of the unstable manifolds of 
periodic orbits are the decisive factors in determining the far-field pattern. Surprisingly, we find 
that strongly chaotic ray dynamics is compatible with highly directional emission in the far-field. 
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I. INTRODUCTION 

Spherical, cylindrical and disk-shaped dielectric cavi- 
ties have been of interest as compact, high-Q optical res- 
onators to be used in micro-lasers and integrated optics 
applications |M l2l| . The high-Q modes of such devices are 
"whispering gallery modes" (WG) which circulate inside 
the boundary confined by total internal reflection. Due 
to their intrinsic rotational symmetry such devices re- 
quire additional symmetry-breaking elements to couple 
in a directional manner and to optimize the Q- values for a 
given application. Some time ago it was shown by Nockel, 
Stone and Chang^J I12I fl3| that smooth deformations 
of such spherically symmetric cavities (termed ARCs for 
asymmetric resonant cavities) still have anisotropic WG 
modes which now have intrinsically directional emission 
and Q-values which were tunable by the degree of de- 
formation introduced. The optical physics of such res- 
onators is non-trivial and interesting because the ray 
dynamics in such a case is partially chaotic and stan- 
dard real-space ray-tracing is not helpful in analyzing 
their properties. Instead, phase space methods taken 
from non-linear classical dynamics such as the surface-of- 
section method (see below) allowed a much clearer pic- 
ture of the physics, leading to qualitative predictions for 
the high intensity emission directions from quadrupole- 
deformed ARCs |fl LL1 UJ. This phase space picture and 
associated ray simulations were shown to agree semi- 
quantitatively with exact numerical calculations of the 
linear resonances of quadrupole ARCs |rj ITU . 

The essence of the Nockel-Stone ray model for ARCs 
was that emission from deformed WG modes could be 
viewed as refractive escape of rays which were initially 
trapped by total internal reflection, but which, due to 
their chaotic dynamics, could diffuse chaotically in an- 
gle of incidence, %, until they reach the critical angle, 
sin% c = 1/n (n is the index of refraction of the res- 
onator, assumed to be surrounded by air) and refracted 
out. One might naively assume that such chaotic ray dy- 



namics would generate a fluctuating and pseudo-random 
emission pattern, but in fact the ray motion follows a 
dominant flow pattern in the phase space favoring escape 
at certain points on the boundary and in certain direc- 
tions in the far-field. Both numerical experiments and 
later experimental measurements on ARC lasers found 
highly directional emission patterns in the far-field. For 
the quadrupole and related ARCs it was argued that the 
flow pattern was approximately describable as rapid mo- 
tion along adiabatic invariant curves (which could be cal- 
culated from knowledge of the boundary shape) and slow 
diffusion in the transverse direction. However major de- 
viations from the flow pattern would occur in the vicinity 
of stable periodic ray orbits for reasons to be discussed 
in detail below. As the shape of the adiabatic curves 
and the location of stable and unstable periodic orbits 
is quite sensitive to boundary shapes, this theory pre- 
dicted a rather dramatic sensitivity of the emission pat- 
terns from ARC micro-lasers to the shape of the bound- 
ary. In the current work we explore a range of interesting 
ARC shapes, both within and outside the range of earlier 
theory. One of our goals is to test the earlier theoreti- 
cal predictions experimentally for the first time with a 
controlled series of boundary shapes for polymer ARC 
microlasers. The shapes we consider are quadrupole, el- 
liptical and hexadecapole ARCs (see definitions, caption 
of Fig. We find that; I) There is a remarkable and re- 
producible difference in the lasing emission patterns from 
ARC lasers with very similar boundary shapes. 2) The 
basic difference between chaotic (quadrupole) and non- 
chaotic (elliptical) ARC emission patterns is in agree- 
ment with the predictions of Nockel and Stone based 
on the adiabatic model. 3) Nonetheless the persistence 
of highly directional emission patterns for highly de- 
formed quadrupolar ARCs is inconsistent with the adi- 
abatic model and is a quite surprising experimental dis- 
covery. 4) A new theoretical model is proposed which 
attributes the high emission directions observed for the 
chaotic shapes to the flow pattern produced by the unsta- 
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ble manifolds of short periodic orbits; this flow pattern 
differs significantly from the adiabatic model. This model 
is shown to explain both the persistence of narrow emis- 
sion peaks in the quadrupole at high deformation and 
the major shift in emission directionality at large defor- 
mation for ARCs with hexadecapole deformation. It also 
predicts that completely chaotic boundary shapes, such 
as the stadium, can nonetheless exhibit highly directional 
emission. 

Further below (section IIV|) we will review the basic 
ideas of mixed phase space, which describes systems with 
a mixture of chaotic and regular dynamics. At that 
point we will explain in detail the adiabatic model of 
Nockel and Stone and how it fails to account for the emis- 
sion data for the highest deformations of resonators with 
mixed phase space. However we review here just a few 
basic concepts to put the current experimental and theo- 
retical work in perspective. At maximal radial deforma- 
tions of order 15% from a reference circle the quadrupole 
and hexadecapole resonators have primarily chaotic ray 
dynamics and this implies that most initial conditions 
corresponding to total internal reflection (TIR) belong to 
ray trajectories which eventually will strike the bound- 
ary below TIR and hence will escape rapidly by refrac- 
tion. The prior model posited a correspondence between 
a set of totally-internally-reflected initial conditions on 
rays and WG modes of the deformed resonator. This 
correspondence, based on adiabatic invariants, will be 
described below. The correspondence is essentially ex- 
act for the ellipse but is only a rough approximation for 
the quadrupole. Once ray initial conditions are chosen, 
the emission patterns can be calculated using straight- 
forward ray simulations as described in^J U| and a 
qualitative understanding of the emission patterns is pos- 
sible based on properties of the phase-space flow. Us- 
ing this picture of phase space flow, Nockel and Stone 
predicted 0, 0,0]; that quadrupole resonators with de- 
formations in the range of 10-12% and index of refraction 
n = 1.5 would emit very differently than elliptical res- 
onators with the same major to minor axis ratio. (Sim- 
ilarly, but less relevant here, they showed that a n = 2 
quadrupole resonator of exactly the same shape would 
emit differently from the n = 1.5 case). Specifically the 
elliptical resonators would emit from the points of high- 
est curvature on the boundary roughly in the tangent 
direction (90°) while the quadrupoles of the same index 
would emit at roughly a 35° —45° angle to the major and 
minor axes. It was argued that the origin of this effect 
is the presence in the quadrupole of a stable four-bounce 
periodic ray orbit which prevents emission from the high- 
est curvature points in the tangent direction, an effect 
they termed "dynamical eclipsing" [Til Il3| . This finding 
was supported by numerical solutions of the linear wave 
equation for the quasi-bound states and their far-field 
emission patterns. Mode selection and non-linear lasing 
processes were not treated in the theory. A subsequent 
experiment on lasing droplets by Chang et al.|6( was suc- 
cessfully interpreted as strong evidence for such dynam- 
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FIG. 1: Cross-sectional shapes of micro-pillar resonators 
studied: a) The quadrupole, defined in polar coordinates 
by R = Ro(l + e cos 20), b) The ellipse, defined by R = 
Ra(l + ((1 + e) 4 - 1) sin 2 0)- 1/2 and c) The Quadrupole- 
Hexadecapole, defined by R = Ro(l + e(cos 2 4> + |cos 4 <jf>)) 
all at a deformation of e = 0.12. Note that all shapes have 
horizontal and vertical reflection symmetry and have been de- 
fined so that the same value of e corresponds to approximately 
the same major to minor axis ratio. \ (Fig. b) is the angle 
of incidence of a ray with respect to the local normal. In a) 
and c) we show short periodic orbits ("diamond, rectangle, 
triangle") relevant to the discussion below. 

ically eclipsed lasing modes. However this experiment 
was less direct than desirable for two reasons. First, the 
droplet was a deformed sphere with many possible lasing 
modes, most of which were not of the two-dimensional 
type considered in references [III fl^| : it was argued, but 
not experimentally shown, that the 2-d "chaotic" modes 
dominated the lasing emission. Second, the droplet shape 
deformation was not controlled and could not be manip- 
ulated to turn the effect on and off. Recently, Lacey 
et al.0 reported an experiment on nearly spherical res- 
onators where they addressed the former but not latter 
problem. The experiments reported here remedy both of 
these shortcomings. First the lasers are deformed cylin- 
ders and the lasing modes are truly two-dimensional. 
Second we have fabricated the boundary shapes using 
a mask to conform to the desired cross-sectional profile. 
Hence we can directly compare e.g. quadrupolar and el- 
liptical ARC lasers and observe the presence or absence 
of the dynamical eclipsing effect over a range of deforma- 
tions. 



II. EXPERIMENTAL DATA 

The experiments we report were performed on differ- 
ently shaped dye (DCM)-doped polymer (PMMA) sam- 
ples that are fabricated on top of a spin-on-glass buffer 
layer coated over a silicon substrate via a sequence of mi- 
crolithography and O2 reactive ionic etching steps. The 
effective index of refraction of these microcavities is 1.49, 
much lower than for earlier experiments performed us- 
ing a similar set-up on GaN, where the index of refrac- 
tion is n = 2.650, [T3 - They are optically pumped by 
a Q-switched Nd:YAG laser at A = 532 nm incident 
normal to the plane of the micropillar. Light emitted 
from the laser is imaged through an aperture subtend- 
ing a 5° angle and lens onto a ICCD camera which is 




FIG. 2: Two-dimensional display of the experimental data showing in false color scale the emission intensity as a function of 
sidewall angle cj> (converted from ICCD images) and of the far-field angle 6 (camera angle) . Columns from left to right represent 
the quadrupole, ellipse and quadrupole-hexadecapole respectively. Insets show the cross-sectional shapes of the pillars in each 
case (for definitions see Fig. 0. The graphs at the bottom show the far-field patterns obtained by integration over <j> for each 
9, normalized to unity in the direction of maximal intensity. The deformations are e = 0.12,0.16,0.18,0.20 (red, blue, black 
and green respectively) 



rotated by an angle 6 in the far-field from the major 
axis. A bandpass filter restricts the imaged light to 
the stimulated emission region of the spectrum. The 
ICCD camera records an image of the intensity profile 
on the sidewall of the pillar as viewed from the angle 
(j> which is converted from pixels to angular position 
<f>. In this paper we studied micro-cavities with ellip- 
tic, quadrupolar and quadrupolar-hexadecapolar shape 



of an average radius R = 100/im (see formulas in Fig. 2] 
caption). Each shape was analyzed at eccentricities of 
e = 0.12, 0.14, 0.16, 0.18 and 0.20. 

In Fig. El we show the experimental results in a color 
scaled ICCD image. The two angles are the sidewall an- 
gle 4> (for the horizontal coordinate) and the camera an- 
gle 9 (for the vertical coordinate). We omit the data for 
e = 0.14 deformation as it indicates no effects not cap- 
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tured by the data at the other deformations. To obtain 
the far-field pattern with respect to 9 we integrate the 
image over all sidewall angles <f>. The boundary image 
field is calculated by integrating over 9 for each tj>. As 
insets we show the exact shape of each of the microcavi- 
ties. Although the shapes appear very similar to the eye, 
we find dramatic differences in the far-field emission pat- 
terns, which in the case of the ellipse vs. the quadrupole, 
persist over a wide range of deformations. Specifically, 
the far-field emission intensity for the quadrupole ex- 
hibits a strong peak at 9 — 34° — 40° which remains 
rather narrow over the observed range of deformations. 
Over the same range of deformation the boundary im- 
age field for the quadrupole changes substantially and 
does not exhibit one localized point of emission. In con- 
trast, the ellipse emits into the 9 ~ 90° direction in the 
far-field, but with a much broader angular intensity dis- 
tribution, while the boundary image field remains well- 
localized around </> ~ 0° (the point of highest curvature 
in the imaged field). For the quadrupole-hexadecapole 
we see a far-field directionality peak which shifts from 
9 ~ 90° to 9 ~ 30° and an almost constant boundary 
image field. Thus we see three qualitatively different be- 
haviors for the three shapes studied over the same range 
of variation of the major to minor axis ratios. 

Several different samples with the same boundary 
shape were measured in each case and confirmed that 
the basic features of this data set just described are re- 
produced within each class (with small fluctuations) ^(| . 
This shows that the effects measured are a property of 
the boundary shape and not of uncontrollable aspects of 
the fabricationprocess. Moreover the theoretical calcula- 
tions, which we will present next, are based on uniform 
dielectric rods with the ideal cross-sectional shape spec- 
ified by the mask; therefore the agreement of these cal- 
culations with the measurements also confirms that the 
differences are due to controllable shape differences. 



III. THEORETICAL CALCULATIONS 

In Figs. 13151 we compare the experimental results for 
the far-field emission patterns for the three shapes mea- 
sured at e = 0.12,0.18 to two theoretical models, one 
based on solutions of the wave equation and the other 
based on simulations of ray escape. The agreement in 
both cases is quite good. We briefly summarize here the 
two models used. 

It is conventional in treating high-Q laser cavities to 
approximate the lasing mode by the appropriate reso- 
nant mode of the passive cavity (the cavity in the ab- 
sence of pumping and gain) 18] . In the current case we 
can model the passive cavity as a 2-d resonator because 
the penetration depth of the pump laser is of the order 
of /jm and amplification due to reflection from the side- 
wall occurs primarily within the horizontal plane of the 
resonator. The wave-equation for the time-harmonic so- 
lutions in 2-d within the uniform dielectric reduces to 




FIG. 3: Far- field intensity for the quadrupole with e = 0.12 
(a) and 0.18 (b). The green curve is the experimental result, 
blue the ray simulation and red a numerical solution of the 
wave equation. The ray simulation was performed starting 
with 6000 random initial conditions above the critical line and 
then propagated into the far-field in the manner described 
in the text. The numerical solutions selected have kRo = 
49.0847 - 0.0379i with a Q = -2Re[kR ]/lm[kR ] = 2593.05 
and kR = 49.5927-0.0679i with Q = -2Re[kR ]/lm[kR } = 
1460.72 for e = 0.12 and 0.18 respectively. 



the Helmholtz equation for the electric field (TM modes) 
and magnetic field (TE modes) and we can solve the 
Helmholtz equation by imposing Sommerfeld boundary 
conditions (no incoming wave from infinity). We use an 
efficient new numerical method developed for these sys- 
tems described in ref. 20] The method yields a large set of 
resonances over the entire range of Q- values; the real part 
of the wavevector giving the resonance frequency and the 
imaginary part giving its width. Experiments have indi- 
cated that mode selection is complicated in these dielec- 
tric resonators and that there is no simple rule relating 
the observed lasing mode to the Q-value of the mode in 
the passive cavity. Due to the low output coupling high 
Q-modes are not necessarily the observed lasing modes 
in the far-field. Thus from the set of calculated reso- 
nances we choose the resonance which coincides well with 
the observed far-field pattern and has a relatively high 
Q = — 2Re[fc]/Im[fc]. We also confirmed that theoreti- 
cal boundary image data|l7j coincides well with the ex- 
perimental results. Moreover in all cases discussed here, 
there were many resonances which gave good agreement 
with the data, indicating the existence of a robust class 
of modes any of which could be the lasing mode. In 
Fig. 13151 we as red line the numerical far-field by calcu- 
lating the asymptotic expansion of our wavefunction in 



FIG. 4: Far-field intensity for the ellipse with s = 0.12 (a) 
and 0.18 (b). Green, blue and red curves are experiment, ray 
simulation and wave solution. The ray simulation was per- 
formed starting with 6000 initial conditions spread over seven 
caustics separated by Asinx = 0.02 below the critical caustic 
(the caustic that just touches the critical line). The numerical 
wave solutions shown correspond to kRo — 49.1787 — 0.0028i 
with Q = -2Re[kR ]/Im[kR ] = 17481.38 and kR = 
49.2491 - O.OllOi with Q = -2Re[kR }/lm[kR ] = 4488.20 
for e = 0.12 and 0.18 respectively. 



the far-field. Numerical limitations prevent us from per- 
forming the calculations at the experimental values of 
kR ~ 1000 but the major features of the emission pat- 
tern are not sensitive to kRo over the range we can study 
numerically. The finding (discussed next) that we can re- 
produce these patterns from ray escape simulations also 
suggests that the wavelength is not a relevant parameter 
for the features we are studying. With a green line we 
show the experimental results. We see that the numerical 
calculations agree with the measured far-field very well. 

As just noted, the experiment is performed well into 
the short- wavelength limit, and we are motivated to 
develop ray-optical models for the emission, following 
references [ill Il3j . In the ray simulation an initial dis- 
tribution of rays is assumed and each ray is given unit 
initial amplitude. At each reflection on the boundary 
the amplitude is reduced according to Fresnel formulas. 
The outgoing amplitude is recorded in the direction de- 
termined by Snell's law and the reflected ray is followed 
until its amplitude falls below 10~ 4 . To compare to the 
experimental data we collected the transmitted rays in 
5° bins. A subtle issue in the calculations is the choice 
of the initial ray distribution. In Figs. [3] and [S] we show 
the far-field distribution for a randomly chosen set of ini- 



FIG. 5: Far-field intensity for the quadrupole-hexadecapole 
with e = 0.12 (a) and 0.18 (b). Green, blue and red curves 
are experiment, ray simulation and wave solution as described 
in the caption to Fig. |3] The numerical wave solutions 
shown correspond to kRo = 50.5761 — 0.0024i with Q = 
-2Re[kR ]/Im[kR ] = 42573.77 and kR = 49.5642 - 0.0092i 
with Q = -2Re[kR ]/lm[kRo] = 10741.93 for e = 0.12 and 
0.18 respectively. 

tial conditions above the critical angle; in the case of the 
ellipse (Fig. we chose initial conditions appropriate to 
its integrable dynamics (see the discussion below). The 
ray model is found to reproduce the main features of the 
data quite well. In section we show that the far-field 
emission pattern for chaotic shapes is insensitive to the 
initial ray distribution over a wide range. Specifically, 
in Fig. ^5] below we compare the far-field patterns for 
different possible initial distributions, confirming the ap- 
proximate independence of the patterns to this choice. 



IV. NOVEL FEATURES OF THE DATA 

In the previous section we showed that we can repro- 
duce the experimental data with two different theoretical 
models. First, by solving the linear wave equation in the 
dielectric and choosing an appropriate resonance; second, 
by ray escape simulations. This gives us confidence that 
the major differences in the experimental emission pat- 
terns are due to the different shapes of the laser cavities. 
The strong sensitivity of the emission patterns to small 
differences in boundary shape is quite striking and is a 
major result of this study. This sensitivity was predicted 
in the earlier work of references H [H EH and was not 
unexpected. However there are major aspects of the ex- 
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perimental data which are quite surprising even in the 
light of the earlier work on ARCs. In particular, the per- 
sistence of highly directional emission in the quadrupolar 
shapes at quite high deformations was not predicted the- 
oretically and was unexpected for reasons we will now 
discuss. In order to understand the unexpected features 
of the data and to develop principles to predict the emis- 
sion patterns for untested boundary shapes we will re- 
view and extend the phase space approach introduced 
originally by Nockel and Stone. 

To explain the basic concepts it is useful to begin by ne- 
glecting the possibility of refraction out of the resonator. 
The problem of ray motion in a resonator with perfectly 
reflecting walls is equivalent to the billiard problem first 
posed by BirkhoffQ and since widely studied in the non- 
linear dynamics community. The billiard problem is that 
of a point mass moving freely in two dimensions within a 
boundary of general shape defined by perfectly (specular) 
reflecting wall. In such a system classical mechanics is 
reduced to a simple geometric construction, whose prop- 
erties nonetheless turn out to be remarkably complex. 
An important tool to study dynamical systems such as 
billiards was devised by Poincare the Poincare Sur- 
face of Section (SOS). The SOS images the trajectories 
in phase space instead of in real space. For even the sim- 
plest non-trivial dynamical systems (two dimensions, one 
degree of freedom) the phase space is four-dimensional, 
however the SOS is a two-dimensional plot which is a 
stroboscopic image of a set of trajectories as they cross 
a certain hyper-plane in phase space. For billiards this 
image is taken at each reflection from the boundary; the 
coordinates of the SOS are the angular position ((f>) and 
the tangential momentum sin \ (scaled to unit total mo- 
mentum) of each trajectory as it reflects from the bound- 
ary (Fig.^). The topology of the resulting plot gives us 
critical information about the phase space flow in the 
particular billiard of interest. 

When the motion is partially chaotic a real-space trace 
of the trajectory is essentially useless, while the SOS can 
reveal useful information. An example of this is shown 
in Fig. EK). Here we show a SOS for the quadrupole and 
the ellipse with a deformation of e = 0.072; we calculated 
and plotted trajectories arising from ~ 30 different ini- 
tial conditions and iterated them for 500 reflections. We 
see that for the chaotic orbit Fig. ||J\) (c) the real space 
picture does not provide any useful information, whereas 
the SOS shows us an underlying structure; specifically 
that such a chaotic orbit nonetheless avoids large regions 
of phase space in this case. This plot illustrates a typical 
structure for a mixed phase space. Such a structure is 
generic for all smooth deformations of a circular billiard 
although the fraction of phase space which is covered 
by chaotic motion increases with deformation. We will 
not attempt to review the general properties of mixed 
dynamical systems which are described by Kolmogorov- 
Arnold-Moser (KAM)@ theory and by other theorems 
due to Poincare, Birkhoff and Lazutkin. However a key to 
understanding a given billiard is to study its periodic or- 



bits. The "diamond" periodic orbit in Fig. ElA) (b) shows 
up in the SOS as four points which are referred to as fixed 
points in the billiard map. For this deformation the di- 
amond orbit is called stable as motion close to the fixed 
points is bounded in the vicinity of the fixed points. Iter- 
ation of a point close to a stable fixed point will result in a 
1-d set of points and is called quasi-periodic orbit as it fills 
a 1-d curve densely. The set of all the points in the neigh- 
borhood of a stable fixed point that are quasi-periodic is 
called an island. Another type of quasi-periodic motion 
is shown in Fig. [SJY) (a) ; in this case we have motion on a 
"KAM" curve which is similar to the WG modes found in 
the circle and is not associated with any stable periodic 
orbit. The final type of motion which occurs in a mixed 
phase space is shown in Fig. EK)(c). Here we show a 
chaotic orbit which fills a two-dimensional region of the 
SOS densely; hence this is a qualitatively different type of 
motion. Within this region nearby initial conditions will 
separate exponentially in time, (until they are separated 
by a distance of order the size of teh chaotic region) . As 
noted above, the existence of these three types of motion 
is generic for smooth deformations of a circle and hence 
for ARC resonators. 

While the behavior of the quadrupole shown in 
Fig. HJ\) is generic there do exist special billiards that 
exhibit the two extremes of dynamical behavior. One 
limiting case is exemplified by the ellipse billiard whose 
SOS is shown in Fig.|Hj3). For the ellipse all orbits lie on 
1-d curves in the SOS and there are no chaotic regions of 
phase space. This type of system is called integrable and 
has regular dynamics because there exists one constant 
of motion for each degree of freedom. For the elliptical 
billiard these constants are the energy and the product 
of the angular momenta with respect to each focus (de- 
generating to the angular momentum in the case of the 
circular billiard) 0, ^| • The ellipse is the only convex de- 
formation of a circular billiard which is integrable ^{|- A 
recent proof of this long standing conjecture was given 
by Amiran lj. At the opposite extreme is the Buni- 
movich stadium billiard (see inset in Fig. 1141) for which 
it is proven that there exist no stable periodic orbits and 
the entire phase space (except sets of measure zero) is 
chaotic. We will study theoretically the emission from 
stadium-shaped resonators in section IV 1111 

Due to its integrability, phase space flow in the el- 
lipse is particularly simple: every initial condition lies 
on one of the invariant curves given by Eq. JQ| below, 
and the trajectory retraces this curve indefinitely (see 
Fig.|Hj3)). Curves which cross the entire SOS correspond 
to real-space motion tangent to a confocal elliptical caus- 
tic Fig.Hj5)(a); curves which do not cross the entire SOS 
represent motion tangent to a hyperbolic caustic in real 
space Fig.|Hf3)(b). 

Phase space flow in mixed systems is much more com- 
plex and is ergodic on each chaotic region. However a 
key property of mixed dynamical systems is that the dif- 
ferent dynamical structures in phase space are disjoint; 
this implies that in two dimensions KAM curves and is- 
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FIG. 6: The Poincare surface of section for the quadrupole 
A) and the ellipse B) with e = 0.072. The schematics A)(a- 
c) on right show three classes of orbits for the quadrupole, 
A) (a) a quasi-periodic orbit on a KAM curve, A)(b) a stable 
period-four orbit, (the 'diamond'), and A)(c) a chaotic orbit. 
Schematic B)(a, b) show the two types of orbits which exist 
in the ellipse, the whispering gallery type, with an elliptical 
caustic B)(a) and B)(b), the bouncing ball type, with a hy- 
perbolic caustic. 
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FIG. 7: Comparison of the Poincare surface of section for 
the quadrupole and the ellipse with e — 0.12 showing mostly 
chaotic behavior in the former case and completely regular 
motion in the latter. The red line denotes sinxc = 1/n, the 
critical value for total internal reflection; rays above that line 
are trapped and those below escape rapidly by refraction. The 
quadrupole still exhibits stable islands at <f> = 0, tt and sin x = 
sin Xc which prevent escape at the points of highest curvature 
in the tangent direction 



lands divide phase space into regions which cannot be 
connected by the chaotic orbits. This puts constraints 
on phase space flow despite the existence of chaos in a 
significant fraction of the phase space. For small defor- 
mations (~ 5%) most of phase space is covered by KAM 
curves, the form of which can be estimated using an adi- 
abatic approximation ^lj- This approximation gives the 
exact result for all deformations in the case of the ellipse; 
it can be written in the following form: 

sin X (0) = yJl + lSP-lWPfae) (1) 

where k is the radius of curvature along the boundary and 
S is a constant. Plotting this equation for different values 
of S,e gives an SOS of the type shown in Fig. |SJ3). For 
the mixed case, exemplified by the quadrupole billiard in 
Fig. EK) , Eq. QJ describes quite accurately the behavior 
for values of sinx near unity, but doesn't work well at 
lower sin% where chaos is more prevalent. 

Up to now the discussion has been based on the ideal 
billiard, where a ray is trapped indefinitely by perfectly 
reflecting walls. Nockel and Stone realized that in this 
short wavelength limit the dielectric resonator would be- 
have very similar to the ideal billiard for angles of in- 
cidence above sinx c = 1/n (i.e. the evanescent leakage 
could be neglected), however for incidence below sin^c 
rays would rapidly escape by refraction according to Fres- 



nel's law. If a correspondence could be made between a 
set of initial conditions for rays and a set of solutions of 
the wave equations (modes) then the emission pattern 
could be calculated by propagating those rays forward 
in time and allowing them to escape into the far-field 
according to SneU's and Fresnel's laws. For the case of 
the ellipse there was an obvious correspondence between 
a set of rays chosen on a given invariant curve of the 
ellipse and a set of solutions. This correspondence can 
be formalized using the eikonal method of Keller 8] . The 
difficulty was in finding the correspondence in the chaotic 
regions of phase space of generic shapes, for which there 
exist no invariant curves. 

In refs.0,^] a model was proposed based on an adi- 
abatic theory to describe the ray-wave correspondence in 
the generic case of mixed dynamics. Initial conditions 
on the adiabatic invariant curves, e.g. of the quadrupole 
were chosen. Due to the presence of chaos in the true 
dynamics, rays initially on such a curve would diffuse in 
phase space until they escaped by refraction. The re- 
sulting emission pattern can be calculated by ray simula- 
tions of the type we have presented above. Moreover this 
model led to qualitative predictions about the emission 
patterns without doing any simulations. The adiabatic 
invariant curves for the quadrupole have their minimum 
values of sin x at the points of highest curvature on the 
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boundary = 0, ±tt, just as they do in the ellipse. If 
the diffusion in phase space is sufficiently slow, emission 
would be near these points of highest curvature and at 
the critical angle, i.e. in the tangent direction, as in the 
ellipse. This reasoning held as long as the escape points 
sinx = 1/n, = 0, ±7r occurred in the chaotic region 
and were reachable from the totally-internally-reflected 
region of sin x > 1/n. However for n — 1.5 and deforma- 
tions around 10%, these points are enclosed by the stable 
island corresponding to the four-bounce "diamond" orbit 
and due to the disjoint nature of the dynamics, "chaotic" 
rays cannot escape there. Instead they will escape at 
higher or lower values of leading to a large change in 
the emission pattern from that of the ellipse with similar 
minor-major axis ratio. This is the phenomenon termed 
"dynamical eclipsing" , and it was predicted to occur for 
the n = 1.5 quadrupole at e ~ 0.12 some time ago [TTl[T3j . 
Our experimental data confirms this prediction for the 
e = 0.12 case. Fig. contrasts the phase space for the 
ellipse and the quadrupole for this deformation. The is- 
land associated with the stable diamond orbit is smaller 
than at e — 0.072, but still present for the quadrupole; 
there is no such island at any deformation for the ellipse. 
Note that in the experimental data for the quadrupole at 
e = 0.12 we do not see a bright spot at the boundary at 
= 0, consistent with the dynamical eclipsing model in 
which the island structure forces the chaotic WG modes 
to emit away from the point of highest curvature. In con- 
trast the bright spot in the ellipse which emits to 9 = 90° 
clearly is a t = for e = 0.12. Thus the adiabatic model 
of refs. 0, El does seem to provide a reasonable descrip- 
tion of the data for e — 0.12 and the observed dramatic 
difference between the ellipse and quadrupole shapes is 
as predicted. 

The earlier work on ARCs did not look extensively at 
deformations above e = 0.12 for the case of low index 
materials such as polymers. The belief was that the adi- 
abatic model would become questionable at higher de- 
formations as the phase space became more chaotic and 
the ray motion departed from the adiabatic curves very 
rapidly. A natural expectation was that due to increased 
chaos the emission patterns in the far-field would become 
less directional and more pseudo-random. More specif- 
ically, for the n — 1.5 quadrupole one finds that the 
"diamond" orbit becomes unstable at e w 0.1369 (the 
associated islands shrink to zero) and one would not ex- 
pect highly directional emission at higher deformations 
(see inset in Fig. |8J. Thus a plausible extrapolation of 
the adiabatic model suggests a steady broadening of the 
quadrupole emission with deformation with at least some 
significant emission in the tangent (90°) direction. Our 
experimental data strongly contradicts this expectation, 
as the observed emission patterns remain peaked around 
35° and do not broaden at all up to e = 0.20. A similar 
analysis would show that the adiabatic theory provides 
no qualitative explanation for the switch in the far-field 
directionality for the quadrupole-hexadecapole at high 
deformations. Thus we are motivated to look for a model 




FIG. 8: Poincare surface of section for the quadrupole with 
e = 0.18. The red line indicates the critical angle of incidence. 
The green diamonds indicate the location of the fixed points 
of the (now) unstable "diamond" orbit and the blue squares 
the fixed points of the unstable rectangular orbit. In the inset 
we show the trace of the monodromy (stability) matrix (see 
Eq. for the diamond orbit versus deformation. When the 
magnitude of the trace of the monodromy matrix is larger 
than two its eigenvalues become real, the periodic motion be- 
comes unstable and the associated islands vanish. For the 
diamond this happens at e — 0.1369 (see dashed vertical line 
in the inset) and the simple dynamical eclipsing picture of 
Fig- El does not apply at larger deformations. 



of the phase space flow which can explain the persistence 
in highly directional emission at high deformations where 
almost all of the phase space is chaotic. This model will 
be presented in the next section. 



V. SHORT-TIME DYNAMICS AND UNSTABLE 
MANIFOLDS 

At high deformations chaotic diffusion is fast and rays 
tend to escape rapidly even if they are initially well con- 
fined (i.e. far away from the critical angle). The adia- 
batic model reviewed predicted at which regions on the 
boundary rays are more likely to escape. How accurate is 
this model as a description of the phase space flow? One 
easy thing to check is whether the short-term flow pre- 
serves the discrete symmetries of the system; the shape 
of the adiabatic curves are determined by the shape of 
the boundary, via its curvature (see Eq.Q and thus must 
preserve all the billiard symmetries (e.g. horizontal and 
vertical reflection axes for the quadrupole). A straight- 
forward ray simulation shown in Fig. [!|] reveals that the 
short-term dynamics typically breaks these symmetries. 
We begin a set of rays on an adiabatic curve, which has 
reflection symmetry around = 0, and after 50 itera- 
tions the SOS density has developed a strong asymmetry 
around the origin. Thus the flow cannot be accurately 
seen as simple diffusion between adiabatic curves. On 
the other hand the resulting SOS density is very struc- 
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tured and does not look like a realization of independent 
random walks. We shall see below that the boundaries of 
the high density regions are set by the unstable manifolds 
of a few short periodic orbits. 
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FIG. 9: Simulation of short-term ray dynamics for the 
quadrupole with deformation e = 0.12, shown in blue. The 
three red curves are examples of adiabatic curves for the 
quadrupole; the blue points are the result of iterating 8,000 
starting rays for 50 bounces starting on the middle of the three 
red curves. The blue points show a complex structure which 
is asymmetric around the line <j> = in contrast to the SOS 
and the adiabatic curves which have reflection symmetry. 



Another indication that the emission patterns we ob- 
serve are not closely linked to adiabatic curves is ob- 
tained by looking at the dependence of the ray emission 
directionality on initial conditions. The original work 
of Nockel and Stone began the ray simulations with an 
ensemble of rays uniformly distributed on the adiabatic 
curves. This is obviously correct for the ellipse as one can 
calculate resonances which cover uniformly the different 
invariant curves. However, as already noted, the ellipse 
is special since chaotic diffusion is absent. In Fig. 1101 we 
compare ray emission patterns arising from three differ- 
ent choices of initial conditions in the quadrupole: uni- 
form on the adiabatic curve (black dashes), initial con- 
ditions localized on the unstable fixed points associated 
with the rectangular periodic orbit shown in Fig. ^(ma- 
genta dashes), and finally, initial conditions chosen ran- 
domly above the critical angle for escape (blue). These 
three quite different choices all lead to similar far-field in- 
tensity patterns, in good agreement with the experimen- 
tal measurements (green). More generally, we found that 
the ray simulations are quite insensitive to the choice of 
initial conditions, as long as a significant fraction of the 
rays are started within the chaotic sea. (Ray bundles 
only started in an island would obviously lead to differ- 
ent results.) We shall propose an explanation for this 
insensitivity in the next section. 




FIG. 10: Ray simulations of the far-field emission patterns 
for the quadrupole with e = 0.12 (a), e = 0.18 (b) with dif- 
ferent types of initial conditions. The blue curve is the result 
of choosing random initial conditions about the critical line 
sin x = l/ n , the black dashed curve is for initial conditions on 
the adiabatic curve with minimum value at the critical line. 
The magenta curve is for initial conditions localized around 
the unstable fixed point of the rectangle periodic orbit. In 
each of the ray simulations 6000 rays were started with unit 
amplitude and the amplitude was reduced according to Fres- 
nel's law upon each reflection, with the refracted amplitude 
"collected" in the far-field. The green curve is the experimen- 
tal result; clearly all three choices give similar results in good 
agreement with experiment. 



VI. UNSTABLE MANIFOLDS 

By definition a non-linear system which is partially 
chaotic, such as our billiards, generates ergodic and mix- 
ing motion on a finite fraction of the phase space for long 
times. It is convenient to analyze the dynamics in the 
SOS and not in continuous time; in this case the dynam- 
ics is described by a discrete map. One can get a good 
idea of the short term dynamics of such a map in selected 
regions of the SOS by linearizing it in the neighborhood 
of unstable fixed points of the billiard map (correspond- 
ing to unstable periodic orbits in real space). If we take 
the initial position and direction/momentum of one ray 
at the boundary to be [s,u) = (jfi, sinx) we define the 
map which projects the ray to the next position and di- 
rection to be 

T: (0, sin X ) - 0', sin *')■ (2) 
A set of fixed points of order N is defined by 

T N (0, sin x) = (0, sin X ). (3) 
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We can propagate an initial ray corresponding to a 
small deviation from the fixed point values by linearizing 
the map around the fixed point. 

(ds'(s,u) du'(s.u) \ 
Osf s ,u) (0 ' Sin 
du du / 

(4) 

The nature of the nearby motion can then be character- 
ized by calculating the eigenvalues and eigenvectors of M . 
For Hamiltonian flows M is always an area-preserving 
map, i.e. detM = 1. The matrix M is also known as the 
monodromy, or stability matrix. The eigenvalues can be 
either complex on the unit circle or purely real and recip- 
rocal to each other. If the eigenvalues are complex, the 
fixed points are stable (elliptic) and nearby points oscil- 
late around the fixed points tracing an ellipse in the SOS. 
In this case the long-time dynamics is determined by the 
linearized map to a good approximation. In the case of 
real eigenvalues there will be one eigenvalue with mod- 
ulus larger than unity (unstable) and one with modulus 
less than unity (stable) and there will be two correspond- 
ing eigendirections (not usually orthogonal). In the sta- 
ble direction, deviations relax exponentially towards the 
fixed points; in the unstable direction deviations grow 
exponentially away from the fixed points. Generic devia- 
tions will have at least some component along the unsta- 
ble directions and will also flow out along the unstable 
direction. Therefore, in a short time generic deviations 
move out of the regime of validity of the linearized map 
and begin to move erratically in the chaotic "sea" . Hence 
the linearized map is not a good tool for predicting long 
time dynamics in a chaotic region of phase space. How- 
ever, in open billiards, rays will escape if they wander 
away from the fixed points into the part of the chaotic 
sea which is below the critical angle for total internal 
reflection. Therefore we find the unstable eigenvectors 
of the short periodic orbits useful in predicting ray es- 
cape. For the short periodic orbits in our shapes it is 
possible to calculate the matrix M giving the linearized 
map around all of the short periodic orbits. Thus we can 
calculate the eigendirections and determine the unstable 
directions analytically. 

Before we can explain the characteristic asymmetric 
pattern found in Fig. we need to introduce the more 
general concept of the non-linear stable and unstable 
manifolds of the fixed points. Consider first the stable 
eigendirection defined by the the stable eigenvector of the 
monodromy matrix for a given fixed point. This is a line 
in the SOS passing through the fixed point such that any 
initial condition on that line will come arbitrarily close to 
the fixed point as t — > oo; similarly the unstable direction 
is a line of initial conditions which approaches the fixed 
point arbitrarily closely as t — > — oo. At a certain dis- 
tance in phase space from the fixed point the linearized 
map no longer holds and a deviation exactly on the line 
corresponding to the unstable eigenvector will not tend to 
the fixed point asymptotically; however nearby this line 
there exist such points which form a generalized curve 



known as the unstable manifold of the fixed point. Each 
unstable fixed point has associated with it stable and un- 
stable manifolds which coincide with the eigendirections 
as one passes through the fixed point. Note that for infe- 
rrable systems there is only one asymptotic manifold for 
both past and future and it coincides exactly with the 
invariant curves, which can be calculated analytically in 
some cases (e.g. the ellipse). For the non-integrable case, 
e.g. the quadrupole, we can only calculate the eigendi- 
rections near the fixed point analytically and must trace 
out the full manifolds numerically. As the unstable man- 
ifolds deviate further from the fixed points, generically 
they begin to have larger and larger oscillations. This is 
necessary to preserve phase space area while at the same 
time have exponential growth of deviations. This tan- 
gling of the unstable manifolds has been used to devise 
a mathematical proof of chaotic motion [Tsj. 

One can argue qualitatively that the unstable mani- 
folds of the short periodic orbits ought to control the ray 
escape dynamics at large deformations. The manifolds of 
short periodic orbits are the least convoluted as they are 
typically the least unstable; hence the unstable direction 
is fairly linear over a large region in the SOS. A typical 
ray will only make small excursions in phase space until it 
approaches one of these manifolds and then it will rather 
rapidly flow along it. If the direction leads across the 
critical line for escape, that crossing point and the por- 
tion just below will be highly favored as escape points 
in phase space. We can check this qualitative argument 
with a few simple numerical experiments. 

In Fig. \llh ) we show the results of a short time itera- 
tion of a uniform random set of initial conditions above 
the critical line in comparison to the unstable mani- 
folds of the various relevant short periodic orbits. Note 
that the different unstable manifolds must fit together 
in a consistent manner and cannot cross one another; if 
they did such a crossing point would define a ray which 
asymptotically in the past approaches two different sets 
of fixed points, which is not possible. Because of this 
non-crossing property the unstable manifolds define just 
a few major flow directions in the SOS. We see clearly in 
the simulation that the actual short-time flow of random 
points in the chaotic sea is controlled and approximately 
bounded by these unstable manifolds. To further support 
our statement that the general motion in phase space is 
governed by the unstable manifolds of these short orbits, 
in Fig. lllb ) we propagated an arbitrary but localized set 
of initial conditions and confirmed that they are stretched 
along and parallel to nearby unstable manifolds. Thus it 
appears that for the highly deformed case the phase space 
flow of a generic ray is much better predicted by simply 
plotting these manifolds. 

As a confirmation that these manifolds do control es- 
cape we perform a further ray simulation for the "open" 
billiard. We propagate, as before, an ensemble of rays 
with a uniform random distribution above the critical 
angle. As we have done in calculating the ray emission 
pattern, we associate to every starting ray in the sur- 
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FIG. 11: (a) Ray simulations of short-term dynamics for 
random initial conditions above the critical line, propagated 
for 10 iterations, plotted on the surface of section for the 
quadrupole with e = 0.18. The areas of the SOS covered are 
delineated very accurately by the unstable manifolds of the 
short periodic orbits which are indicated in the schematics 
at right. These manifolds are overlaid in the figure with ap- 
propriate color coding, (b) Flow of phase space volume in 
the surface of section of the quadrupole with e = 0.18. A 
localized but arbitrary cloud of initial conditions (red) is it- 
erated six times to illustrate the flow. The initial volume is 
the circle at the far left, successive iterations are increasingly 
stretched by the chaotic map. The stretching clearly follows 
closely the unstable manifold of the rectangle orbit which we 
have plotted in blue. 



face of section an amplitude which decreases as the ray 
propagates forward in time according to Fresnel's law (if 
the point falls below the critical line) . Instead of follow- 
ing the refracted amplitude into the far-field, in this case 
we plot the emitted amplitude onto the surface of sec- 
tion, as shown in Fig. 112b ). The emission amplitude is 
almost completely comfined within the two downwards 
"fingers" created by the unstable manifold of the four- 
bounce rectangular orbit. As noted earlier, the availabity 
of the two-dimensional data obtained from the imaging 
technique (see Fig. gives us a unique ability to re- 
construct the emitting part of the lasing mode both in 
real space and momentum space directly from experi- 
mental data. It is therefore possible to check directly 
this ray simulation in phase space against experimental 
data. The intensity data is sorted into intensity pixels 
according to both its sidewall location (the angle cj> from 
which emitted intensity originated) and its far-field an- 
gle, which by geometric considerations and SnelPs law 
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FIG. 12: a) Ray simulation of emission: emitted ray ampli- 
tude (color scale) overlaid on the surface of section for the 
quadrupole with e = 0.18. b) Far-field intensity from exper- 
imental image data Fig. [5] projected in false color scale onto 
the surface of section for the quadrupole with e = 0.18. The 
blue line is the unstable manifold of the periodic rectangle 
orbit. In green we have the line of constant 34° far-field (see 
the discussion in section IV 111 . Absence of projected intensity 
near <j> — ±7r in (b) is due to collection of experimental data 
only in the first quadrant. 

can be converted to the internal angle of incidence sin%- 
Therefore we can project this data "back" onto the SOS 
for emission. In Fig. 112b ) we show this projection for 
the same deformation as in Fig. 112b ): we find remarkable 
agreement between the projected data and the ray sim- 
ulation. We note that this is a much more demanding 
test of agreement between theory and experiment than 
simply reproducing the experimental far- field patterns. 



VII. RAY DYNAMICAL EXPLANATION OF 
THE EXPERIMENTAL DATA 

In the last section we established that typical rays 
above the critical angle escape by following closely the 
unstable manifolds of the short periodic orbits. This 
leads to a ray escape probability which is relatively local- 
ized in the surface of section (Fig. 112b)). However despite 
the non-random character of this escape, there is still a 
significant spread of angles of incidence for escape. In 
fact the spread of escape angle we see in Fig. IT21 would 
lead to an angular spread of nearly 80° in the far-field 
if all the escape occurred from the same point on the 
boundary. However as we see from Fig. 1121 the point of 
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FIG. 13: Ray emission amplitude (color scale) overlaid on 
the surface of section for the quadrupole-hexadecapole with 
e = 0.12 (a) and e = 0.20 (b). Solid blue and magenta curves 
are the unstable manifold of the diamond orbit (a) and of the 
unstable rectangular orbit (b) . In green and turquoise we plot 
the line of constant emission in the 75° and 30° directions in 
the far-field. 



escape and the angle of incidence are correlated and vary 
together according to the shape of the unstable mani- 
fold. Because the boundary is curved, different angles of 
incidence can lead to the same angular direction in the 
far-field. It is straightforward to calculate the curves of 
constant far-field for a given shape; for the quadruple at 
e = 0.18 and for the peak observed emission angle of 34° 
this curve is plotted in green in Fig. 1121 The curve tends 
to lie remarkably close to the unstable manifold. There- 
fore we find that the curvature of the boundary tends to 
compensate almost completely for the dispersion in the 
angle of incidence at escape. 

As a further test of the explanatory power of plot- 
ting the unstable manifolds, we can use this method to 
explain the large shift in the far-field directionality in 
the quadrupole-hexadecapole (QHD) with deformation 
shown in Fig. [21 The QHD shape is an interesting con- 
trast to the quadrupole as the diamond and rectangle 
periodic orbits interchange their roles. For the QHD at 
small deformations the diamond orbit is unstable and 
gives no islands whereas the rectangle is stable and gives 
rather large islands at e — 0.12 (see Fig. 113b')'). Since 
there is no island at the point of highest curvature for 
this deformation the original adiabatic theory would pre- 
dict emission from the point of highest curvature approx- 
imately in the tangent (90°) direction. We see in Fig.0 
that we indeed have such behavior experimentally. The 
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FIG. 14: Far-field emission patterns for the stadium with 
e = 0.12 (a) and e = 0.18 (b). The blue curve is the 
ray simulation and the red a numerical solution of the wave 
equation; no experimental data was taken for this shape. 
The ray simulation was performed with random initial con- 
ditions exactly as in Fig. |3] The numerical solutions were 
for resonances with kR = 50.5401 — 0.043H with Q — 
-2Re[kR]/Im[kR] = 2342.71 and kR = 48.7988 - 0.1192i 
with Q = -2Re[kR]/lm[kR] = 818.83 for e = 0.12 and 0.18 
respectively. The inset shows the shape of the stadium; it is 
defined (inset (b)) by two half circles with radius one and a 
straight line segment of length 2e. 



same prediction would come from looking at the unstable 
manifolds; in Fig. 113b ) we find the maximum ray escape 
amplitude comes from near cf> — and is bounded by 
the unstable manifold of the diamond orbit. The rele- 
vant unstable manifolds rearrange as the deformation is 
increased to e — 0.20. The stable rectangle bifurcates 
at e ~ 0.1115 into two stable, period four, "parallelo- 
gram" orbits. As the deformation increases the islands 
associated with these orbits move closer to <j> = 0, and al- 
though the islands themselves are quite small, they cause 
the unstable manifolds of the diamond to become steeper 
(as the fixed point of the diamond orbit sits right between 
the two period four islands) . Due to this the other mani- 
folds also become steeper and at e — 0.20 Fig. 113b) shows 
that the unstable manifold of the rectangular orbit dic- 
tates the flow of the escaping rays, which now are emitted 
into the 9 = 30° direction in the far-field for essentially 
the same reason as in the quadrupole. Thus the QHD 
is a shape which behaves like the ellipse at low deforma- 
tions and as the quadrupole at high deformations; this 
can be attributed to the evolution in the geometry of its 
unstable manifolds. 
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FIG. 15: Ray emission amplitude (color scale) overlaid on 
the surface of section for the stadium with e = 0.12 (a) and 
e = 0.18 (b). Solid blue curve is the unstable manifold of the 
periodic rectangle orbit. The green curve is the line of con- 
stant 55° (a) and 48° (b) emission direction into the far-field. 
The thick black lines mark the end of the circle segments of 
the boundary and coincide with discontinuities in the mani- 
folds. 



VIII. DIRECTIONAL EMISSION FROM 
COMPLETELY CHAOTIC RESONATORS 

The existence of highly directional emission for the 
highly deformed quadrupole (e = 0.20) suggests that 
the slow diffusion in phase space, characteristic of mixed 
systems, is not essential to get this effect. Therefore 
we decided to study theoretically resonators for which 
the corresponding billiard is completely chaotic and for 
which there exist no stable periodic orbits at all. The 
Bunimovich Stadium (see inset in Fig. I14|) . mentioned 
above, was a natural choice due to its similarity to the 
quadrupole. As before we did both ray escape simula- 
tions and numerical solutions of the wave equation. In 
Fig. we show our predictions. We find again highly di- 
rectional emission with a peak direction (~ 55°) slightly 
shifted from the quadrupole; the narrowness of the far- 
field peak in the stadium is comparable to that of the 
far-field peak in the quadrupole. We can associate this 



peak with the slope and position of the manifold of the 
unstable rectangular orbit in the stadium, Fig. ll5h ,). The 
noticable shift between the e = 0.12 and e — 0.18 defor- 
mation (see inset in Fig. I14|l originates from the change 
in the slope of the unstable manifold of the rectangu- 
lar orbit, Fig. 115b). The discontinuities of slope in the 
unstable manifolds of the periodic orbits in the stadium 
result from its non-smooth boundary. These results indi- 
cate clearly that a fully chaotic dielectric resonator can 
nonetheless sustain highly directional lasing modes. It 
would be interesting to test this in future experiments. 



IX. SUMMARY AND CONCLUSIONS 

We have reported experimental data from polymer 
micro-pillar lasers with different deformations of the 
cross-section from circular symmetry. The far-field emis- 
sion patterns were anisotropic and in most cases highly 
directional. The anisotropy was stable and reproducible 
and was dramatically different for different boundary 
shapes of similar major to minor axis ratio. These dif- 
ferences were reproduced by ray and wave solutions of 
the ideal passive cavity. The differences in the emission 
patterns were explained by reference to the different dy- 
namics of rays trapped within each resonator. The pos- 
sibility of highly directional emission from quadrupole 
resonators with partially chaotic ray dynamics was pre- 
dicted earlier, but was found to be much more robust 
in the experiment than previously suspected. The ear- 
lier adiabatic model^lj could not explain this robustness 
and a new model based on the geometry of unstable man- 
ifolds was proposed and tested. It was able to explain the 
current data and by extension predicts that fully chaotic 
resonator shapes such as the stadium would also exhibit 
highly directional emission. 
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